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Abstract 


It  is  necessary  to  develop  detailed  numerical  forecasting 
techniques  that  accurately  specify  the  process  of  wave  generation 
because  the  need  for  accurate  wave  forecasts  is  increasing.  The 
growth  of  the  wave  spectrum  can  be  approximated  by  the  combined 
effect  of  two  mechanisms  to  a point  where  the  crest  of  the  wave 
breaks  on  a nonlinear  effect  becomes  dominant. 

One  mechanism  is  the  Phillips  resonance  mechanism  for 
wave  generation,  and  the  other  is  the  exponential  growth  mechanism, 
originally  proposed  by  Miles.  Recently,  this  second  growth  mech- 
anism was  revised  by  Phillips  in  terms  of  the  wave-induced  atmo- 
spheric perturbations.  Calculations  suggest  that  this  new  mechanism 
explains  the  available  observations,  which  had  previously  been 
unexplainable  from  the  original  instability  growth  mechanism. 

To  verify  the  theory,  various  spectra  estimated  from 
wave  observations  are  used,  along  with  measurements  provided  by 
fetch  limited  field  experiments  made  by  Snyder  and  Cox.  After 
some  assumptions  about  energy  dissipation  are  applied,  the 
spectral-growth  equation,  which  is  a function  of  frequency,  wind 

speed,  and  duration  (or  fetch),  is  obtained.^ 

\ 

The  partially  developed  sea-spectral  shapes  provided 
by  the  spectral  growth  equation  for  several  different  wind  speeds 
are  shown  and  compared  to  previous  work.  The  proposed  spectral 
growth  shows  faster  development  in  an  early  growing  stage  (or 


short  fetch),  and  slower  growth  in  the  last  stages.  It  takes  from 
35  to  40  hours  or  700  to  1000  nautical  miles  to  reach  a fully- 
developed  sea  for  a 40  knot  wind.  The  growth  rate  of  the  signi- 
ficant wave  height  is  about  the  same  as  the  results  of  Sverdrup 
and  Munk. 

The  proposed  expression  for  the  spectral  growth  is 
applied  in  a numerical  forecasting  computer  program  to  the  North 
Atlantic  for  the  month  of  December  1959.  The  time  history  of  the 
hindcasted  significant  wave  height  and  the  spectral  density  for  low 
frequencies  are  compared  to  the  British  weather  ships'  observed 
data.  The  wave  height  shows  good  agreement  with  the  observed, 
with  the  bias  calculated  to  be  less  than  +2  feet  and  the  RMS 
error  less  than  5 feet.  Excellent  verification  is  also  obtained 
with  Argus  Island  data. 
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1 . Introduction 

In  the  modern  high-speed  computer  era  of  today,  the  techniques 
for  forecasting  ocean  waves  are  being  studied  not  only  for  their  purely 
scientific  content  but  also  for  the  vast  wealth  of  information  that  can 
be  supplied  to  such  practical  problems  as  shipping,  ship  designing, 
and  forecasting  natural  disasters.  The  past  several  years  have  shown 
an  ever-increasing  interest  in  all  phases  of  oceanographic  activity. 

This,  then,  necessitates  a greater  understanding  of  the  natural  pheno- 
mena called  oceanic  surface  waves.  Consequently,  the  study  and 
forecasting  of  ocean  waves  must  not  be  limited  to  only  one  ocean,  but 
to  all  three  of  the  world  oceans.  The  use  of  large,  high-speed  com- 
puters has  made  rapid  and  precise  determinations  of  oceanic  wave 
characteristics  feasible. 

The  actual  study  of  ocean  wave  forecasting  can  be  said  to 

have  begun  during  World  War  JI  by  Sverdrup  and  Munk  (1947).  This 

study  introduced  the  concepts  of  significant  wave  height,  H,  wave 

duration,  t,  fetch,  F,  steepness  (wave  height/wave  length)  6,  and 

wave  age  (wave  phase  speed/wind  speed)  (3,  as  well  as  relationships 

2 

between  non-dimensional  quantities  such  as  gF/u  , gt/u,  6 and  (3 

(where  g = gravitational  acceleration,  and  u = wind  speed).  If  either 
2 

gF/u  or  gt/u  were  determined,  then  6 and  (3  can  be  found.  6 
and  (3  were  then  used  to  determine  the  wave  height  H,  and  the 
wave  speed  C.  Several  years  later,  Bretschneider  revised  this 
method. 

In  1952,  the  successful  application  of  wave  spectra  concepts 
to  ocean  waves  was  accomplished  in  studies  of  randomness  of  the  sea 


2 


surface  (Pierson  and  Marks,  1952).  This  allowed  the  idea  of  a com- 
posite ocean  wave  consisting  of  different  frequencies  and  amplitudes 
to  be  incorporated  into  a forecasting  scheme.  Prior  to  this,  the 
classical  method  of  studying  ocean  waves  considered  each  wave  to  be 
a sinusoid.  Also  about  this  time,  Longuet-Higgins  (1952)  made  a 
statistical  study  of  ocean  wave  properties  and  presented  statistical 
relationships  for  the  significant  wave  height,  the  average  height  of 
the  one -tenth  highest  waves  and  the  mean  wave  height  in  terms  of  the 
variance  of  a wave  record.  In  the  following  years,  Neumann  (1953) 
proposed  a spectral  form  and  co-cumulative  spectra  (CCS-curves)  that 
were  derived  from  data  taken  visually  on  board  a ship  and  off  Long  Branch, 
New  Jersey. 

By  incorporating  these  ideas,  Pierson,  Neumann  and  James  (1955) 
were  able  to  develop  a new  wave  forecasting  technique  based  on  spectral 
concepts.  This  new  forecasting  technique  was  very  different  from 
the  significant  height-significant  period  method  in  that  a great  amount 
of  information  could  be  gained  about  the  nature  of  ocean  waves. 

Thus  in  the  middle  1950's,  two  schools  of  thought  existed  for 
the  "technique"  of  forecasting  ocean  waves--the  "significant  height" 
school  and  the  "spectral"  school. 

With  the  improvement  of  wave  recording  instrumentation 
(i.  e.  , Tucker,  1956),  more  accurate  observations  became  available. 

These  data  analyses  were  then  followed  by  many  proposals  of  the 
growth  and  spectral  form  of  ocean  waves  (i.  e.  , Darbyshire,  1955, 

1959;  Bretschneider,  1959).  Several  attempts  were  made  to  measure 


t- 


the  directional  spectrum.  Cote  et  al  (1960),  for  instance,  worked  with 


3 


data  obtained  from  pairs  of  aerial  stereophotographs  while  Longuet- 
Higgins  et  al  (1963)  dealt  with  data  obtained  from  a buoy.  But  most 
of  the  studies  were  directed  towards  finding  the  form  of  the  one- 
dimensional spectrum.  With  the  host  of  proposals,  came  many 
differences  of  opinion.  For  instance,  Darbyshire  concluded  that  the 
significant  wave  height  was  proportional  to  the  square  of  the  wind 
speed  whereas  Neumann  (1953)  concluded  that  the  proportionality 
was  to  the  2.5th  power.  These  discrepancies  were  discussed  at  the 
Conference  on  Ocean  Wave  Spectra  (Ocean  Wave  Spectra,  1963).  A 
study  of  these  discrepancies  was  made  by  Pierson  (1964)  in  terms  of 
the  variation  of  wind  speed  with  height.  The  use  of  a wind  profile 
to  normalize  the  wind  speeds  to  the  same  height  resulted  in  a better 
agreement  among  the  proposed  relations  between  significant  wave 
height  and  wind  speed. 

While  a wave  is  growing  due  to  wind  effects,  the  high  frequency 

range  of  the  spectrum  grows  more  rapidly  than  the  low  frequency 

band  because  of  greater  energy  input.  At  the  same  time,  dissipation 

caused  by  wave-wave  interaction,  breaking  waves,  and  viscous  effects 

is  also  increasing  and  thus  a so-called  ’’equilibrium  range"  of  the 

spectrum  is  rapidly  approached.  Phillips  (1958)  found  the  general 

character  of  this  equilibrium  range  and  proposed  an  expression  of 
2 5 

the  form  ag  /u>  (where  a = constant,  g = gravitational  acceleration, 
and  a)  = angular  frequency)  for  the  spectrum  in  this  range.  Pierson 
and  Moskowitz  (1964)  then  proposed  a spectral  form  for  a fully 
developed  sea  that  was  based  on  the  similarity  theory  of  Kitaigorodskii 
(1961)  and  data  taken  by  British  weather  ships. 


The  development  of  high-speed  electronic  computers  allowed 
the  field  of  ocean  wave  forecasting  to  make  significant  advances. 
Large-scale  numerical  wave  forecasting  techniques  were  attempted 
(Baer,  1962).  The  generation  and  growth  processes  of  waves  which 
were  not  understood  to  any  appreciable  extent  at  this  time  were  in- 
corporated by  use  of  energy  growth  tables.  It  was  assumed  that  the 
transient  spectrum  had  the  same  shape  as  the  fully  developed  sea 
spectrum  for  lower  wind  speeds  that  produced  the  same  significant 
wave  heights.  Up  until  this  time,  forecasting  was  done  manually 
by  first  employing  the  various  graphs. 

The  demand  for  oetter  prediction  accuracy  led  to  the  study  of 
wave  generation,  growth  and  the  partially  developed  sea  spectrum. 
The  study  of  wave  generation  and  growth  dates  back  to  the  Helmholtz 
instability  theory  (Lamb,  1932).  This  Kelvin -Helmholtz  model  could 
not  adequately  explain  wave  growth  since  it  required  stronger  winds 
than  are  actually  observed  to  generate  waves. 

Jeffreys  (1925),  in  his  studies  on  the  initial  phases  of  wave 
formation,  introduced  the  concept  of  a sheltering  coefficient.  This 
was  an  expression  for  the  increase  of  pressure  on  the  windward 
side  and  decrease  of  pressure  on  the  leeward  side  when  wind  passed 
over  a pre-existent  wave.  The  pressure  difference  distribution  was 
then  used  to  determine  a number  called  the  sheltering  coefficient 
empirically.  Experiments  on  a solid  wavy  model  later  showed 
Jeffreys'  results  to  be  about  an  order  of  magnitude  too  large. 
Sverdrup  and  Munk  (1947)  suggested  that  a smaller  value  could  be 
obtained,  and  they  took  into  account  the  tangential  stress  acting  on 


the  sea  surface  in  addition  to  the  normal  stress.  This,  however,  leads 
to  a contradiction  since  a tangential  stress  causes  waves  to  be 
rotational  and  actual  observations  have  shown  waves  to  be  almost 
ir  rotational. 

It  should  be  noted  that  the  sheltering  coefficient  concept  assumes 
pre-existent  waves.  If  the  sea  surface  were  calm  with  no  wavy 
motions,  then  waves  could  not  grow. 

A coupling  mechanism  between  turbulent  components  of  the  at- 
mosphere and  sea  for  wave  generation  from  a glassy,  calm  sea  was 
later  proposed  by  Phillips  (1957).  This,  in  essence,  said  that  a 
resonance  between  the  air-sea  system  could  occur  when  a component 
of  the  surface  pressure  distribution  moved  at  the  same  speed  as  a 
free  surface  wave  of  the  same  wave  number.  This  resonance  would 
then  cause  waves  to  be  generated  even  from  calm  sea  surface  con- 
ditions. 

In  another  attempt  to  overcome  the  sheltering  coefficient  dif- 
ficulties, Miles  (1957)  put  forth  another  theory  after  computing  the 
amplitude  and  lag  of  the  pressure  fluctuations  over  waves  in  terms 
of  properties  of  the  logarithmic  wind  profile.  This  theory  said 
that  the  mean  rate  of  energy  transferred  from  the  parallel  shear 
flow  to  the  surface  wave  is  proportional  to  the  curvature  of  the  wind 
profile,  and  inversely  proportional  to  the  slope  of  the  wind  profile 
at  the  height  where  the  mean  wind  velocity  is  the  same  as  the  phase 
speed  vf  the  wave  component  under  investigation. 

In  further  work  in  this  field,  Miles  (1959)  obtained  a solution 
to  the  Orr-Sommerfeld  equation  after  including  the  effects  of  a 
viscous  term.  A combined  resonance-instability  model  for  wave 


generation  was  also  worked  out  by  Miles  (I960).  Lighthill  (1962) 
later  gave  a physical  interpretation  to  this  model. 

Snyder  and  Cox  (1966)  made  a series  of  field  experiments 
(based  on  these  new  generation  theories)  to  study  the  growth  of  a 
wave  component  in  a limited  fetch  area.  Their  results  confirmed 
the  effects  of  the  two  proposed  mechanisms  although  a difference  in 
the  magnitude  of  the  growth  rate  was  found. 

With  the  existence  of  these  two  new  theories,  the  study  of  wave 
forecasting  can  now  be  carried  out  from  a new  point  of  view,  namely 
that  of  wave  generation.  An  attempt  to  employ  these  two  theories 
in  terms  of  spectral  concepts  has  already  been  made  (Inoue,  1966). 
Large-scale  forecasting  schemes  are  also  attempting  to  incorporate 
these  ideas  (Pierson,  Tick,  and  Baer,  1966). 

Thus,  more  precise  and  detailed  forecasting  of  deep  water  ocean 
waves  can  now  be  accomplished  by  utilizing  wave  generation  concepts. 
The  purpose  of  this  study,  therefore,  is  to  predict  the  character 
and  spectrum  of  deep  water  ocean  waves  at  a given  location  and  time 
by  use  oi  a computer-based  technique  utilizing  the  concepts  of  wave 
generation  mentioned  above.  . 

2.  Generation  and  growth  of  the  waves 

In  order  to  forecast  ocean  waves,  a study  must  be  made  to 
determine  how  spectra  of  the  waves  can  grow.  When  a component 
of  the  spectrum  is  considered,  that  spectral  component  can  be  ex- 
pressed in  the  linear  form  in  the  early  stage,  for  nonlinear  effects 
at  this  stage  may  be  neglected. 
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■—  S(f,  t,  3T)  = A(f,  u(t,  x ))  + B(f,  u(t,  x ))S(f,t,  X ) 


(2.  1) 


where  S(f,  t,  x)  is  the  spectral  density  and  A(f,  u(t,  x ) ) and 

B(f,  u(t,xT))  are  unknown  quantities  to  be  determined.  If  these 
two  functions,  A and  B,  are  known,  the  growth  of  the  spectra 
can  be  computed  for  an  early  growing  stage.  In  the  old  model 
(Inoue,  1966),  these  two  quantities  were  estimated  as 


A(u)  = 1.4  x 10'8u3 


B(f,  u)  = 6.27(u/c)2e"°  °17{c/u)  f 


(2.2) 


where  u is  the  wind  speed  (knots)  at  19.5  m anemometer  height 
and  c is  the  phase  speed  of  the  component  f.  The  A term  was 
just  approximated,  and  it  must  be  a function  of  frequency.  The  B 
term  is  also  changed  in  terms  of  the  new  theory. 

The  pressure  fluctuations  on  the  sea  surface  can  be 
decomposed  into  two  categories.  Cne  represents  the  contribution 
of  the  atmospheric  turbulent  eddiss  and  the  other  is  the  contribution 
induced  by  the  air  flow  over  the  sea  surface.  The  sea  surface 

can  be  expressed  by  the  form,  d£  (k  ) e , where  k 

is  a wave  number  and  d£(k  ) is  a random  complex  quantity. 

The  air  over  that  sea  surface  induces  a pressure  distribution 

given  by  dw(k,  t)e^*  C ^ . If  the  phase  of 

d^(k,  t)  is  the  invariant  with  respect  to  d£(k  ),  the  mean 


(dw)  is  a pressure  field  induced  by  a sea  whose  wave  number  k 
is  the  same. 

The  relation  between  the  pressure  distribution  and  the  sea 
surface  wave  amplitude  was  shown  (Miles,  1957,  1959)  in  the  idealized 


form: 


p = Re  {(a  + ib)  pauJ  k£} 


where  a and  b are  real,  non-dimensional  quantities  depending  on  the 
wind  profile,  pa  is  the  density  of  air,  u^  is  u#./kQ  (u^  is  the  friction 
velocity  and  kQ  is  von  Karman's  constant),  and  £,  is  the  sea  surface. 

Phillips  used  a more  realistic  description  which  is  ex- 
pressed by  a random  wave  amplitude  d£(k  ) and  an  induced  random 
pressure  dw(k,  t). 


<dw)  = (v  + ip.)  p c kd£(k  ) 


where  v and  p.  are  non-dimensional  coupling  coefficients,  pw  is 
the  water  density,  c is  the  phase  speed,  and  the  notation,  ( ) , 

represents  an  averaging  over  time  in  a frame  of  reference  moving 
in  the  direction  k at  the  appropriate  phase  speed. 

Thus  the  pressure  field  is  a combination  of  two  parts:  one 
is  induced  from  the  sea  configuration  which  is  expressed  as  above; 
and  the  other  is  a random  quantity  that  comes  from  the  atmospheric 
turbulence.  This  turbulent  pressure  supplies  energy  over  a wide 
frequency  region,  and  an  induced  pressure  caused  by  the  sea  surface 
configuration  feeds  back  to  the  sea  and  also  contributes  to  the  growth 


of  the  wave. 

Phillips  (1966)  solved  the  equation  governing  the  growth  of  a 


wind  generated  sea  under  the  assumption  that  the  water  is  irrota- 
tional  and  that  a linealized  free  surface  boundary  condition  could  be 
used.  The  surface  displacement  spectrum  ^(k  , t)  was  obtained.  When 
time  t is  larger  than  an  integral  time  scale,  and  wt  >>  1,  ^(k  , t) 

is  expressed  by  using  the  wave  number,  frequency  atmospheric 
pressure  spectrum  l!(k  ,cj)  as 


<t(iT , t)  = Amhtii 


2 2 

p c 
rw 


P O) 


(2.5) 


Miles  (I960)  also  obtained  this  relationship.  If  poo  t is  much  less 
than  one,  then  sinhpwt  — poot  so  that  (2.5)  becomes 


V(k  , t) , sM-kjA  t 

0 c 
• w 


(2.6) 


Equation  (2.6)  shows  that  the  wave  spectrum  grows  linearly 
with  time,  t,  and  the  atmospheric,  turbulent  fluctuation  spectrum. 
This  atmospheric  turbulence  spectrum  shows  its  maximum  at  the 
frequency  co  = k • u^  , where  uc  is  the  convection  velocity.  There- 
fore the  sea  surface  spectral  component  becomes  dominant  at  the 
frequency,  f,  which  is  the  same  frequency  as  the  maximum  turbu- 
lent spectral  component.  This  is  a typical  resonance  phenomenon. 

If  poot  is  much  larger  than  one,  then  sinhpoot  — e^^/Z  , so 
that  equation  (2.  5)  can  be  simplified  to 


*{k  . t)  = Tr-?i-k-2’^ 

2pwc  pto 


(2.  7) 


This  equation  (2.7)  shows  that  the  wave  spectrum  grows  exponentially 


with  time  t in  the  growing  stage  (root  >>  1,  and  this  is  an  instability- 
growth  mechanism. 

Thus,  there  are  two  different  stages  in  the  generating  and  grow- 
ing sea  spectrum.  First  the  sea  starts  to  grow  from  the  glassy,  calm 
sea  in  terms  of  the  resonance  mechanism,  and  then,  later  the  growth 
by  an  instability  mechanism  is  involved.  The  latter  mechanism  domi- 
nates the  former,  and  the  sea  grows  exponentially.  Jeffreys  tried  to 
determine  the  atmospheric  pressure  distribution  induced  by  this  wave 
motion  empirically  by  introducing  a sheltering  coefficient  concept. 

These  mechanisms  can  be  interpreted  in  equation  (2.1)  as  follows. 
First  the  sea  is  calm,  and  there  is  no  spectral  energy  over  the  entire 
frequency  range,  and  no  contribution  is  made  by  the  second  term  B S. 
Then  the  sea  grows,  and  the  second  term  plays  an  important  role.  In 
equation  (2.1),  the  first  term  is  caused  by  the  atmospheric  turbulent 
pressure  tluctuation,  and  the  second  term  is  due  to  the  instability 
mechanism.  This  is  also  evident  when  equation  (2.7)  is  differentiated 
with  respect  to  time  t.  Then, 

^|*(r,t)  = ^(r.t)^  (2.8) 

This  is  exactly  the  same  expression  as  (2.1)  if  the  first  term  A 
on  the  right  hand  side  is  negligibly  small. 

3.  Resonance  mechanism  growth;  the  A term 

The  atmospheric  turbulent  fluctuation  spectrum  is  needed  to  get 
the  surface  displacement  spectrum  Skf  k ,t).  Although  this  study  is 
concerned  with  one -dimensional  spectra,  it  will  start  with  the  three- 
dimensional  pressure  spectrum,  because  the  data  used  were  obtained 
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with  a narrow  band  sensor.  Let  p(x,y,  t)  be  the  pressure  at  the 
point  (x,  y)  at  time  t;  then  the  correlation  function  R(§»  t|,  t)  is 

R(£>  rj,  t)  = (p(x,  y,  t)  p(x+£,  y + r)»  t + T)>  (3.1) 

where  £ and  t)  are  deviations  from  x and  y along  their  axes, 

and  t is  a time  lag.  Then  the  pressure  spectrum  can  be  shown  to  be 

II(k*co)=— / d | dr  R(T,T)cos(iT-  T- wt)  (3.2) 

( 2t r)  - oo 

where  £ is  the  vector  form  for  £ and  q . 

Priestly  (1966)  made  a field  measurement  of  the  longitudinal, 
lateral  and  diagonal  correlation  of  the  atmospheric  pressure  over 
grass,  and  this  result  is  used. 

The  space -time  covariance  R(a,  p,  t)  is  transformed  in  a 
Fourier  manner,  and  the  spectral  density  is 

. oo 

K(|,T),«)-  ^ /R(S,-n,T)e‘lwT  dr  (3.3) 

-OO 

Priestly  assumed  that  the  spectral  density  can  be  expressed 
as 

K(e,T),oo)  ^w)e“Yiei"6|Tll_iKe 

where  $(w),  y>  and  6 are  determined  from  the  experiment,  and 
K = w/uc  . uc  is  a convection  velocity.  Then  equation  (3.2) 


becomes 


n(iT,w)  =Re  /— i-y  / d2  £[$(&)  e-^1  "6lT1|]el(k  * ^~K^1 

lC2tr)  -oo  J 


Then 


n(ir,u)=ReJ^  "d6«ieil-''+i(k^-K)!  j d^e|ri|l-6+itkyi3} 
L(2tt)  -OO  -00  J 


where  and  k^  are  wave  numbers  of  the  x and  y components.  This 
can  be  expressed  in  the  following  form: 


n(k  , 


w)=i£y  ‘ 
(2*)  L 


(2tt)2  L y^  + (k  cos0  - K)  J L 6 2 + k2s 


2 

in  9. 


where  k is  defined  as  w = gk  and  0 is  an  angle  from  the  wind  direction. 


According  to  Priestly's  experiment,  y and  6 were  determined  as 
y = 0.33K1*28 

n on  (3-6) 


6 = 0.52K 


Priestly's  observation  shows  some  tendency  for  y to  be  constant 
for  K < 0.02  (m  ^),  but  in  this  study  y is  treated  as  a function  of  K 
as  defined  in  equation  (3.6). 

<p(u)  was  found  to  be  expressed  as 

<P  (w)  = $4,(w)  (3 

where  $ is  a scaling  factor  and  »]/  (w)  is  determined  from  Priestly's 


measurement,  and 


+(«)  — 


In  order  to  determine  the  scaling  factor  $ , the  field  experi- 
ment data  obtained  by  Snyder  and  Cox  (1966)  were  used,  and  $ was 
assumed  to  be  proportional  to  the  fourth  power  of  the  wind  speed, 
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$ oc  const*  u . 

The  one -dimensional  spectrum  is  then  obtained  by  integrating 
(3.5)  with  respect  to  the  angle  9 from  -tt/2  to  tt/2  , and  the  A 
term  is  written  in  the  form 


7 , 3 t/2  _ 

A(f,  u)  = ■ -<%-  / n(kjW)de 

p g -Tr/2 
6 


(3.9) 


After  putting  all  the  constants  together  as  A*  , equation  (3.9)  becomes 


ir/2 

A(f,u)=  / 


*.  ,5.25  .2.25 


-■Jr/2  fl 


2 , , 2.5  2 


d0 


(3.10) 


^0  +<ksine)  +(k  cose-J  ] 


where  w = 2irf,  u is  the  wind  speed  at  a certain  anemometer  height, 

2 ❖ 

k = o>  / g , and  A will  be  determined  later  from  observational 
data. 


4.  Instability  mechanism  growth;  the  BS  term 

Miles  (1957,  1959)  studied  this  instability  type  growth  in  terms 
of  inviscid  laminar  flow.  Phillips  (1966)  took  into  account  the  un- 
dulatory  turbulent  flow  due  to  the  waves.  This  atmospheric  per- 
turbation due  to  the  sea  is  quite  different  from  the  other  mechanisms 
and  affects  the  wave  growth  factor  p strongly.  The  Phillips  theory 
is  summarized  _=  an  important  part  of  this  study  in  this  section. 

The  wind  velocity  in  a moving  frame  with  speed  c is 
expressed  by 


u = U(z)  + u'(x,  y,  z,  t)  + U (x,  z)  - c 

P 


for  a y-direction  averaged  wave  surface  ( £ ) = a cos  kx  , where 


(4.1) 


U(z)  is  the  mean  wind  velocity,  u'(x,  y,  z,t)  is  u random  fluctuation 
from  the  mean,  and  U^(x,  z)  is  a perturbation  induced  by  the  sea 
configuration. 

Following  the  derivations  of  Phillips,  assume  the  atmosphere 
to  be  an  incompressible  fluid;  then  the  stream  function  ijj  is  defined 
as 

3i|< /3z  = U cos  9 + Up  - c 
30i/3x  = -W 

P 

Because  and  W are  periodic  along  the  x-coordinate,  the  stream 
line  of  the  mean  motion  is 

Z ikx 

* = / {u(0  cos  8 - c}  d£.  + 4*(z)  e = const. 

z 

m 

where  is  the  height  of  the  matched  layer  which  is  U(z^) cos0  = c. 

Then  investigate  how  the  momentum  flux  shifts  from  the  atmo- 
sphere to  the  wave.  Use  the  equations  of  motion  and  average  over 
the  y-direction.  The  mean  stress  balance,  then,  is  obtained  by 
taking  the  complete  average  on  the  derived  equation.  That  is. 
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is  the  sum  of  the  wave-induced  stress,  t , and  the  turbulent  stress, 
r^_  . Only  the  first  term  remains  in  inviscid  laminar  flow.  How- 
ever, Phillips  showed  that  this  is  also  a sufficient  approximation 
even  if  the  flow  is  turbulent. 

Then,  considering  the  wave  induced  Reynolds  stress  by  using 
the  product  where  q is  the  vorticity  of  the  wave-induced 

perturbation  and  W is  the  vertical  component  of  the  induced  per- 
turbation. 


fiW  = (~U  - ~ W )W  = (U  W ) + (U  2-  V/  2) 
p 9z  p 9x  p p 9z  p p 2 9x  p p 


which  is  derived  under  the  assumption  of  incompressibility, 

9U  /9x  = - 9W  /9z  . 

P P 

Integrate  (4.6)  with  respect  to  z,  and  take  an  x-average. 

Then 


p f nW  dz  = - p U W = tw(z) 


If  z = 0,  equation  (4.  7)  becomes 


tw(0)  = P / nw  dz 

O V 


(4-  8) 


The  variations  in  the  mean  vorticity  are  supposed  to  be 


proportional  to  the  mean  vorticity  gradient  and  the  magnitude  of 


the  undulations 


|£2i  cc  | U"(z)  6 cos9! 


(4.  9) 


because  the  basic  flow  has  a mean  vorticitv  distribution  U'(z)cos9. 
The  quantity,  6,  is  the  displacement  of  the  mean  streamline  about 


pro  
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the  average  height,  z,  and  is  derived  from  equation  (4.3). 
Outside  the  matched  layer,  6q  is  expressed  by 

W(z  ) 

6 = = 

° k[U(zj)  cos8  - c] 


(4. 10) 


Therefore, 


r*wp  = m |- 


•U"{z)Wp(z)  cosG  ■ 
U(z)  cos0  - c 


(4.11) 


In  an  inviscid,  laminar  flow,  variations  in  ft  and  W are  exactly 
out  of  phase,  so  that  this  term  must  be  zero.  However,  Phillips  con- 
sidered that  there  is  no  reason  to  suppose  ft  and  W to  be 
uncorrelated  in  turbulent  flow.  So  this  term  cannot  be  neglected 
while  we  neglect  it  in  a laminar  flow  model.  Thus,  outside  the 
matched  layer,  equation  (4.10)  survives  and  M is  determined  as 
M =*  1.6  x 10~  2 . 

Inside  the  matched  layer,  the  displacement  of  the  mean  stream 
lines  6m  is  expressed  by 


6 = 
m 


= {w 


4W  (z  ) -,1/2 

a ‘ T*y* » I 


[z  ) cos  9 
^ m 


(4.12) 


The  increment  to  the  wave -induced  Reynolds  stress,  t ,across 

wm 

the  matched  layer  is  derived  from 


. . f-U"  ,.7  2 1 
T = M p vrprr  W }■ 
mw  mr  ^kU1  p Jz_i 


(4.13) 


According  to  Miles'  quasi-lamina r analysis,  = ir  in  the 

above  equation.  This  is  the  form  originally  proposed  by  Miles  (1957). 
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Therefore,  the  total  wave-induced  Reynolds  stress  is  a sum  of 
these  two  equations  in  turbulent  flow. 

, 772 

r (o)  = M pf—iw2!  + Mp 
w'  ' rrv  | kU1  p f ^ 

L v J z=z 

m 

Two  different  values  are  used  for  outside  the  matched  layer  and 
for  inside  the  matched  layer,  respectively,  as  the  approximations  to 
W as  it  appears  in  equation  (4.14). 

Outside  the  matched  layer: 

W = - i N ka  (U(z)  cos0  - c}  e'kzelkx  (4.15) 

Inside  the  matched  layer: 

WP  =-WWe  / {U(Z)cose-c}2e'kzdz  (4.16) 

r ' in'  z 

m 

where  N is  a constant  and  has  a different  value  depending  on 
the  layer  under  consideration,  and  a is  the  amplitude  of  the  free 
surface  wave. 

By  using  the  coupling  parameters  given  in  (2.4),  the  induced 

surface  pressure  distribution  can  be  expressed  as 
2 Ilex 

(v  + iu)p  c kae  . The  imaginary  part,  p,  plays  an  important  role 

W 

in  the  growth  of  the  wave,  and  can  be  expressed  by 

2tw(°) 

(J-  = TT~1  (4*17) 

p c k a 
*w 


,oo  i-uy  w cos  y 

r o £ d2 


ii 


k|  U(z)  cos0  - c | 


'} 


(4.14) 


where  Tw(°)  = ( p (d^/8x),  and  4 *s  the  free  surface.  The 


jjps ■ I V*.*- J" 


combination  of  (4.14),  (4.15),  (4.16),  and  (4.17)  yields 
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pa  , ,-M._N2k2 

d 


”ZT 


w c 


, (-M  IN  K > T,..\  , <x>  _ . w 

^{itv(^L  (^[u“s9-c!  e"  Zdz) 


m m 


2 -2kz  1 

+ M / N (-Uu)dos0  | U cos0  - c|  e dz  l 

Q ■> 


(4.18) 


To  relate  the  curvature  U"  and  the  slope  U1  of  the  wind 
profile  to  the  wind  profile,  we  assume  neutral  atmospheric 
stability.  So  the  wind  profile  becomes  logarithmic  and  is  expressed 
by 


u<*>  - die  = 5T 1”: 


•Z 

o m 


where  u is  the  friction  velocity,  and  k is  a von  Karman  con- 
stant.  From  equation  (4.19),  the  roughness  parameter  zq  is 
expressed  by 


-k  c 
o 

z = z exp-i 

o m r 1 u^cos  0 


Fhillips  used  a relationship  for  zq  suggested  by  Charnock 
(1955)  to  connect  the  roughness  parameter  and  friction  velocity. 
For  aerodynamically  rough  flows,  Charnock's  form  is 


u.„ 

z ex  0.078  — — 
c g 


To  evaluate  equation  (4.18),  N must  be  estimated. 

1 2 

zr  for  z >z  and  N = 1 for  z < z are  used. 

3 m m 

The  value  of  the  growth  rate  p (or  B(f,  u)/f  in  (2.1))  is 


(4.19) 


(4.20) 


(4.21) 
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shown  in  figure  3 by  a dashed  line. 

Kitaigorodskii  and  Volkov  (1965)  gave  data  and  analyses  that 
tended  tc  show  that  the  roughness  parameter  is  not  simply  a function 
of  the  friction  velocity,  but  is  also  a function  of  the  sea  state.  If 
this  is  correct,  the  expression  for  zq  shown  above  is  not  appropriate. 
If  this  type  of  expression  is  allowed,  Kitaigorodskii  and  Volkov  gave 


a different  value  for  the  constant  as 


z — 0.035  — - 
° g 


(4.22) 


As  mentioned  above,  it  is  not  appropriate  to  formulate  a 
relationship  between  the  roughness  parameter  and  the  friction  velocity 
over  the  ocean.  Their  data  also  scattered  widely,  but  a mean  square 
distribution  of  In  zq  fitted  very  well  with  the  Kolmogorov's  test  of 
goodness  of  fit  to  a normal  distribution  in  a given  u^  interval.  If 
this  expression  of  zQ  is  used  in  the  first  term  of  (4.18),  the  curve 
of  the  growth  rate  p is  raised  as  shown  in  figure  3. 

In  the  region  of  u+/c  > 0.06,  the  phase  speed  is  small  and  the 
matched  layer  lies  at  a high  wind  profile  curvature  U",  and  the 
Miles  formula,  the  first  term  in  equation  (4.18),  gives  a large  con- 
tribution to  the  waves.  When  u^/c  decreases,  so  does  the  wind 

ri* 

profile  curvature  and  that  term  becomes  less  dominant.  But  the 
second  term  of  (4.18)  which  is  due  to  the  undulatory  turbulent 
flow  becomes  dominant  and  the  growth  rate  p increases  again. 

This  is  a very  interesting  phenomenon,  and  some  observational 
data  also  show  this  tendency.  The  discussion  of  the  observational 
data  is  given  in  the  next  section. 


5.  Data  and  formulation  for  growth 


To  verify  these  growth  mechanisms,  five  different  kinds  of  ob- 
servations were  used.  The  weather  conditions  were  studied  for  the 
period  from  April  1955  to  March  I960,  and  the  spectra  for  various 
synoptic  situations  were  estimated  from  the  wave  records  obtained  by 
the  British  weather  ships,  "Weather  Explorer"  and  "Weather  Reporter" 
(Moskowitz,  Pierson,  and  Mehr,  1962,  1963,  1965).  These  data  were 
taken  with  the  Tucker  shipborne  wave  recorder  (Tucker,  1956)  at  the 
weather  stations  A,  I,  J,  and  K.  The  location  of  the  ships  are  as 
follows.  Station  A:  62°N,  33°W;  Station  I:  59°N,  19°W;  Station  J: 
52.5°N,  20°W;  Station  K:  45°N,  16°W.  Moskowitz  et  al  picked  460 
wave  records  from  that  period.  The  wave  observations  were  made  at 
three  or  six  hour  intervals,  and  were  approximately  fifteen  minutes 
long.  The  spectra  have  a frequency  range  of  0 ~ 0.333  cycles  per 
second.  The  wind  speed  was  measured  at  19.5  m above  sea  level. 

Wave  records  obtained  at  Argus  Island  were  also  studied.  The 
data  used  were  taken  with  a resistance  wire  wave  recorder  on  Argus 
Island  tower  (32°N,  65°W)  for  the  period  from  20  November  to 
30  November  1961  (Pickett,  1962;  DeLeonibus,  1962).  This  tower 
is  located  25  miles  southwest  from  Bermuda  Island,  and  the  water  is 
approximately  200  feet  deep.  The  observations  were  made  every 
three  hours  and  each  record  was  20  minutes  long.  The  frequency 
covers  0 ~0. 266  cycles  per  second.  The  anemometer  height  was 
42.5  m above  the  sea. 

The  data  obtained  on  FLIP  (Floating  Instrument  Platform) 
were  also  checked.  The  position  of  FLIP  was  approximately 
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39°20'N,  148o20'E  in  the  North  Pacific,  and  the  observations 
were  three  hours  long  and  were  taken  twice  a day.  Sea  surface  wave 
spectra  were  estimated  from  the  records  of  two  vibrotron  trans- 
ducers at  31  m and  88  m below  the  sea  surface  (Snodgrass  et  al, 
1966).  The  data  were  spectrally  analyzed  for  six  30-minute 
pieces  from  three  hours  of  continuous  records  for  the  period  of 
September  1963  (Inoue,  1967). 

These  three  sets  of  data  include  fully  developed  seas,  growing 
seas,  and  decaying  seas.  Only  the  data  from  the  growing  sea,  how- 
ever, are  useful  in  this  study,  and  little  change  of  meteorological 
conditions  between  two  observations  is  preferable.  From  this  point 
of  view,  three  criteria  were  set  at  first:  (1)  significant  wave  height 
had  to  increase;  (2)  the  wind  speed  had  to  change  less  than  two  knots 
within  three  hours;  (3)  the  wind  direction  had  to  change  less  than  45° 
within  the  three  hours  between  two  successive  observations.  When 
the  data  passed  these  criteria,  the  spectra  were  checked  again  as 
follows.  The  spectra  that  contained  dominant  spectral  components 
in  the  low  frequency  swell,  or  two  or  more  big  peaks  were  rejected. 
Under  these  conditions,  only  ten  pairs  of  spectra  were  chosen.  Be- 
sides these  criteria,  one  more  limitation  was  imposed  for  the  pur- 
pose of  this  study.  This  is  (4),  the  spectral  component  of  the 
second  observation  had  to  be  less  than  one-third  of  the  spectral 
component  of  the  fully  developed  sea  for  a particular  frequency. 

The  amount  of  available  spectral  component  growth  data  was  then 
very  limited,  but  these  spectra  and  spectral  components  seem 
to  be  really  growing  stage  spectra. 
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Fetch  limited  wind  wave  studies  were  made  by  Snyder  and 
Cox  (1966)  and  by  Barnett  and  Wilkerson  (1966).  Snyder  and  Cox 
made  an  observation  to  see  how  a particular  spectral  component 
grew  with  fetch  in  the  Bahamas  in  the  spring  of  1963.  Barnett  and 
Wilkerson  took  observations  by  using  an  airborne  altimeter  on  20 
February  1965.  The  flight  was  taken  along  the  direction  SSW 
up  to  190  nautical  miles  off  the  New  Jersey  coast.  These  two  field 
experiments  were  made  mainly  for  the  purpose  of  the  investigation 
of  the  Miles -Phillips  mechanism. 

Snyder  and  Cox  towed  rafts  arranged  17  m apart  downwind 
with  a group  velocity  of  waves  with  a 17  m wavelength,  and  measured 
the  development  of  that  component  wave.  The  wind  speed  ranged 
from  about  10  to  20  knots,  and  the  anemometer  height  was  6.1  m 
above  the  sea  surface.  The  data  were  obtained  for  a total  of  51  hours. 

*{c 

To  evaluate  the  constant  A in  the  A term  (3.10),  the  value 
obtained  by  Snyder  and  Cox  was  used.  Their  wave  recorder  had  an 
effective  beam  width  of  20  degrees.  By  taking  this  width  into 
account  , the  final  form  of  their  result  is  expressed  by 
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where  u is  a wind  speed  in  m/sec  at  19.5  anemometer  height.  (The 

wind  speed  was  corrected  to  this  anemometer  height  from  the  wind  observec 

2 

at  the  lower  anemometer  height),  k = w /9.8. 

The  result  of  this  ‘orm  is  plotted  in  figure  1 which  also  has 
the  value  of  the  A term  estimated  by  Snyder  and  Cox.  The  one- 


dimensional  representation  for  this  form  is  shown  in  figure  2 as  a 

function  of  wind  ■'peed  in  knots.  As  pointed  out  in  Section  3,  y 
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follows  a tendency  toward  being  constant  and  not  a K ’ law  for 
K < 0.02  (m**^).  Then  the  forward  leading  edge  of  the  curve  in 
figure  1 maybe  smaller  and  steeper  for  high  wind  speeds. 

Clearly,  from  the  physical  meaning  of  the  resonance  mechanism, 
the  sea  is  excited  by  an  atmospheric  turbulent  pressure,  and  the 
growth  of  the  A term  is  dominant  at  the  component  that  has  the  same 
velocity  as  the  turbulent  convection  velocity.  The  location  of  the  peak 
on  each  curve  is  about  at  the  frequency  whose  phase  velocity  is  the 
same  as  the  convection  velocity.  In  the  previous  study  (Lnoue,  1966), 
this  A term  was  simply  treated  as  a function  of  wind  speed  over  the 
entire  frequency  range  which  was  not  sufficient. 

In  order  to  investigate  the  growth  with  respect  to  an  instability 
mechanism,  five  different  kinds  of  data  were  checked  as  mentioned 
previously.  Unfortunately,  few  spectra  of  FLIP  for  the  analyzed 
period  passed  the  criteria  described  previously.  The  reason  was  that 
there  were  mostly  moderate  seas  in  September,  and  when  a strong 
wind  was  blowing,  criteria  (2)  and  (4)  were  not  satisfied. 

It  is  difficult  to  detect  the  A term  growth  from  the  spectra  esti- 
mated by  the  British  weather  ship  data,  and  the  Argus  Island  tower 
data.  Therefore,  the  A term  was  considered  to  be  given  by  equation 
(5.1),  and  this  growth  quantity  was  subtracted  from  the  observation 
to  study  the  BS  term  growth.  If  the  value  of  the  A term  is  smaller 
than  the  actual  spectral  component,  the  BS  quantity  can  be  studied 
by  investigating 
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Figure  2 

Resonance  mechanism  growth  A(f,u) 
for  different  wind  speeds. 
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where  So  is  an  initial  spectrum  and  is  a spectrum  after  At  hours. 

In  Section  4,  the  friction  velocity  u_^  was  used  to  evaluate  the 
roughness  parameter  zq.  In  the  actual  ocean,  it  is  a very  difficult 
problem  to  determine  friction  velocity,  and  sometimes  an  approxi- 
mation, u = 9u  ,/k  , has  been  used,  where  k is  von  Karman's 
vo  o 

constant. 

However,  Kitaigorodskii  and  Volkov  (1965)  showed  that  there 
was  a close  relationship  between  the  roughness  parameter  and  the 
sea  state.  When  the  roughness  parameter  is  plotted  against  the 
wave  height,  the  observed  data  scatter  widely,  and  there  is  no  strong 
correlation  between  two  quantities.  However,  the  following  form 
was  found  by  considering  that  the  roughness  parameter  was  dependent 
on  both  the  average  phase  speed  and  the  wave  height. 


zq  = 0,1 20  H exp 
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where  H is  an  average  wave  height,  and  c is  an  average  phase 

speed  of  the  sea.  If  the  relation  of  H = 0.625  H^  is  introduced 
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where  H,  is  the  significant  wave  height,  then 
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z = 0.075  Hi  exp  J - l 
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When  the  atmospheric  condition  is  neutral,  the  wind  profile 
becomes  logarithmic.  The  average  wind  speed  at  height  z is 
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expressed  by  using  this  roughness  parameter. 
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If  the  spectrum  of  that  sea  state  is  known,  an  average  phase 
speed  c can  be  expressed  by 


? = = / ssm 
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where  T is  the  average  period  of  the  sea. 

As  the  above  equation  for  the  friction  velocity  (5.6)  shows, 
u^  is  strongly  affected  by  the  average  speed  c . While  the  waves 
are  growing  at  a certain  wind  speed,  the  spectrum  develops  from 
the  high  frequency  region,  and  then  u is  large  because  of 

/N/  <V 

small  c.  As  the  sea  grows  further,  c becomes  large,  and  u„, 
decreases.  Physically,  this  can  be  interpreted  as  the  fact  that  the 
movement  of  the  wave  system  is  approaching  the  movement  of  the 
atmospheric  system. 

The  friction  velocities  for  the  data  from  the  British  weather 
ships,  Argus  Island  tower,  and  the  airborne  altimeter  were  com- 
puted by  using  the  estimated  spectra  and  equation  (5.6)  for  uv.  The 
spectra  used  were  modified  by  applying  the  ag^/u>^  law  for  the 
equilibrium  range  for  the  high  frequency  region  if  the  available 
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spectra  did  not  cover  the  full  frequency  range.  The  experiment  of 
Snyder  and  Cox  (1966)  was  made  on  a single  component  and  could  not 
reproduce  spectra  correctly.  Therefore,  the  friction  velocity  for 
that  case  was  estimated  from  the  observed  wind  speed  and  the  drag 
coefficient  at  10  m,  Cjq  , proposed  by  Deacon  and  Webb  (The  Sea, 
1962). 

The  observed  results  were  plotted  against  u Jc  to  study  the  BS 
term  growth  in  figure  3.  The  value  on  the  ordinate  shows  B/f 
(or  growth  factor  p).  The  data  obtained  by  Snyder  and  Cox  (1966), 
and  Barnett  and  Wilkerson  (1966)  are  shown  by  + and  • symbols, 
respectively.  The  other  coded  symbols  plotted  are  from  the  British 
weather  ship  data  (Moskowitz,  Pierson,  and  Mehr,  1962,  1965, 

1965)  and  Argus  Island  tower  data  (Pickett,  1963).  Careful  attention 
must  be  given  to  the  data  from  the  Argus  Island  tower  because  the 
geographical  conditions  (the  distance  from  both  the  east  coast  of  the 
United  States  and  Bermuda,  and  the  water  depth)  and  the  small 
difference  between  u and  c might  introduce  some  error.  The 
frequency  appears  to  have  been  high  enough  to  permit  its  gene- 
ration to  full  development  for  the  fetch  and  duration  possible  at 
Argus  Island.  The  data  for  these  observations  such  as  coded 
symbol,  record  number,  wind  speed,  significant  wave  height, 
and  location  are  listed  in  Appendix  I. 

In  figure  3 there  is  a difference  between  the  observed  re- 
sults shown  by  the  coded  symbols  and  the  theoretical  value 
shown  by  the  dashed  and  dotted  lines  [both  lines  are  theoretical 
but  different  roughness  parameters  are  used]  in  the  range  of  high 
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and  low  u^/'c.  In  the  low  u /c  region,  all  the  observed  results 

decrease  as  u /c  decreases,  and  there  is  no  result  close  to  the 
v 

theoretical  values  for  u^/c  < 0.025.  In  the  range  where  uj,./c>0.06, 
the  theoretical  values  are  lower  than  the  observed  results,  and  the 
difference  is  about  a factor  of  3.  When  the  roughness  parameter 
(4. 22) given  by  Kitaigorodskii  and  Voiko/  (1965)  is  used,  the  theo~ 
retical  value  increases,  but  not  by  very  much.  The  difference  be- 
tween the  theoretical  value  and  the  proposed  growth  factor  was  esti- 
mated as  a factor  of  8 in  previous  work  (Inoue,  1966). 

Shemdin  and  Hsu  (1966)  carried  out  an  experiment  in  a water 
tank  on  the  atmospheric  pressure  distribution  in  the  vicinity  of  pro- 
gressive water  waves.  In  the  experiment,  the  growth  of  mechanically 
generated  waves  was  also  studied  and  compared  with  the  theoretical 
result  (Miles’  theory).  Some  of  the  experimental  results  show 
growths  three  times  greater  than  the  theoretical  results,  and  con- 
firm the  observed  results  shown  in  figure  3. 

The  wind  profile  was  treated  as  a logarithmic  profile,  but 
the  profile  depends  on  the  atmospheric  stability,  as  is  well  known. 
Most  data  used  in  this  study  show  a slightly  negative  value  for  the 
air-sea  temperature  difference.  The  air-sea  temperature  differences 
for  'be  British  weather  ship  data  are  listed  in  Appendix  II.  In  Snyder 
and  Cox’s  experiment,  the  conditions  were  nearly  neutral,  and  the 
largest  difference  was  -0.7°C.  Because  the  data  of  Barnett  and 
Wilkerson  were  taken  from  an  airplane,  the  measurement  of 
temperature  w-as  not  made  along  the  flight  track.  However,  it 
was  reported  that  the  difference  was  -6JC~  -10°C  on  the  light  ships 


stationed  near  the  track. 

When  the  atmosphere  is  unstable,  the  friction  velocity  will  be 
higher  tc  an  extent  determined  by  the  Richardson  number.  As  mentioned 
earlier,  the  observed  air-sea  temperature  difference  was  slightly  nega- 
tive, and  the  observed  results  were  plotted  for  neutral  conditions  by 
using  equation  (5.6).  When  the  air-sea  temperature  difference  is  taken 
into  account,  the  plotted  points  of  the  result  might  be  shifted  a litcle 
toward  the  right.  However,  both  sets  of  observations ---the  British 
weather  ship  data  and  the  airborne  altimeter  data--were  taken  during 
a strong  wind,  and  a high  wind  makes  the  Richardson  number 


small.  Therefore,  even  though  the  observed  air-sea  temperature 
difference  was  slightly  negative,  the  growth  rate  was  plotted  against 
the  neutral  conditions  u^  (5.6). 

The  most  important  difference  in  the  new  instability  growth 
mechanism  proposed  by  Phillips  (1966)  is  that  the  mechanism  is  ex- 
pressed by  two  terms--one  is  the  term  originally  proposed  by  Miles 
and  the  other  is  a term  contributed  from  outside  the  matched  layer 
due  to  the  undulatory  turbulent  flow  over  the  waves.  The  latter 
mechanism  makes  the  growth  rate  p (or  B/f)  large  in  the  low  u./c 
region.  However,  the  observed  growth  rate  does  not  show  any  in- 
creasing tendency  as  the  proposed  theoretical  value  increases 
(figure  3). 

One  of  the  most  dominant  results  in  this  study  is  the  "bump" 
existing  in  the  range  of  u^/c  = 0.025  ~ 0.040  in  figure  3.  When  the 
growth  rate  in  the  old  model  (Inoue,  1966)  was  studied,  the  same  kind 


wm 


of  figure  was  plotted  in  order  to  study  the  Miles  mechanism.  And 
the  observed  growth  rate  B/f  also  showed  this  rise  (at  that  time, 
the  quantity,  B/f,  was  plotted  as  a function  of  u/c).  This  rise  and 
fall  toward  decreasing  u^/c  is  not  small.  It  is  quite  significant  be- 
cause the  abscissa  has  a logarithmic  scale.  The  explanation  of  why 
this  hillock  occurs  could  not  be  given  before.  At  first,  the  method  of 
data  selection  was  thought  improper,  but  the  data  were  picked  in  the 
manner  of  the  criteria  mentioned  in  Section  5,  and  the  wind  speeds 
were  also  nearly  constant  during  the  two  successive  observations. 

The  growth  rate  by  the  newly  proposed  process  fits  well  to  the 
observed  data  in  the  range  of  u , /c  = 0.03  ~ 0.05  in  figure  3.  Thus, 
this  "bump1'  can  be  considered  to  be  explained  as  the  result  of  the 
undulatory  turbulent  flow  theory  proposed  by  Phillips  in  that  range 


of  u^/c. 

When  the  friction  velocity  u _ is  constant,  the  spectral  compon- 
ent whose  wavelength  or  phase  speed  is  large  falls  in  the  range  of 
small  u^/c.  The  matched  layer  of  such  a wave  component  is  high 
in  the  atmosphere,  and  the  effect  due  to  the  undulatory  turbulent  flow 
of  that  longer  wave  may  be  considered  to  be  small. 

In  order  to  apply  these  wave  generation  mechanisms,  the  ob- 
served growth  rates  were  used.  The  curve  of  the  growth  factor  B/f 
is  obtained  from  the  plotted  observed  data,  and  is  shown  by  the  solid 
line  in  figure  3.  The  variation  of  B(f,  u^  is  expressed  by 


B(f,u^ -|o.OOi?9e'700C^K/c^"0'031^  4-0.725(_^)  e"°-0004^c/u^2j.  f (5.9) 


The  unit  is  (sec-*). 


As  a form  for  this  type  of  growth,  Snyder  and  Cox  (1966)  and 
also  Barnett  (1967)  estimated  that  it  was  linear  with  respect  to  u/c, 
and  was  zero  for  u/c  < 0.8  ~ 0.9. 

In  equation  (5.9),  the  friction  velocity  formula  proposed  by 
Kitaigorodskii  and  Volkov  is  used.  According  to  their  theory,  u.,. 
is  a function  of  the  sea  state,  namely,  the  average  wave  height  and 
average  wave  speed.  Therefore  the  B term  growth  is  also  a function 
of  the  sea  state.  Or,  stated  another  way,  the  input  energy  for  this 
growth  depends  on  the  sea  state.  When  this  growth  rate  is  only  a 
function  of  u/c,  the  energy  input  of  this  mechanism  is  constant  for 
the  same  wind  speed  for  a particular  frequency.  However,  if  B/f 
is  a function  of  u ^ , the  growth  rate  has  different  values  even  for 
the  same  wind  soeed  and  same  frequency.  Excluding  the  "bump" 
in  the  range  0.03  ~ 0.05  u^/c,  when  the  seas  are  low,  u^/c  is  high, 
and  more  energy  is  fed  in,  generally  speaking.  The  more  the  sea 
grows,  and  the  smaller  u^/c  becomes,  the  less  the  energy  that  can 
be  fed  in.  This  is  just  a general  description,  and  the  "bump"  must 
be  taken  into  account.  However,  this  is  a behavior  different  from 
that  of  the  previous  mode’ 

6.  Spectral  growth 

If  a nonlinear  effect  is  not  involved  in  the  growing  stage, 
equation  (2.1)  can  be  used  to  find  the  generation  and  growth  of  the 
wave  spectra.  As  the  wave  grows,  however,  nonlinear  effects  such 
as  breaking  waves  and  wave -wave  interactions  are  introduced. 

Ocean  wave  spectra  have  been  studied  during  the  last  decade, 
and  many  spectral  forms  and  growth  rates  were  proposed.  Known 
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and  unknown  properties  were  also  discussr 


, . , Neumann  and 


Pierson,  Ocean  Wave  Spectra,  1963).  Ar-,*<ng  those  properties, 
Phillips  (1958)  investigated  breaking  wa\  3s  by  using  similarity  con- 
siderations. As  the  wind  waves  grow,  the  wavy  surface  becomes  un- 
stable at  some  stage.  Then  the  waves  break,  lose  energy,  and  keep 
their  stability.  There  is  some  region  called  the  "equilibrium  range" 
which  attains  a statistical  equilibrium.  In  this  range,  Phillips 
proposed  that  the  spectral  form  can  be  written  as 
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The  energy  transfer  between  different  wave  components  in  a random 
wave  field  was  also  studied  by  Hasselmann  (1962,  1963a,  1963b). 

Equation  (2.  1)  is  linear,  and  the  spectral  component  grows  to 
an  infinite  value  with  time.  But  once  waves  generate,  energy  dissi- 
pation Will  follow  more  or  less  immediately.  Dissipation  has  not 
been  studied  very  much  so  far,  and  little  is  known  about  it. 

There  is,  however,  some  general  agreement  on  the  spectral  growth 
sequence.  The  spectral  component  increases  slowly  after  starting 
from  a zero  initial  condition,  but  dissipation  is  not  dominant.  The 
growth  is  linear  in  the  early  stage;  then  it  becomes  exponential.  This 
is  obvious  from  the  first  and  second  terms  in  equation  (2.1).  When 
this  spectral  component  comes  close  to  the  steady  state,  dissipation 
also  increases,  and  the  growth  is  slowed  down.  Finally,  that 
spectral  component  reaches  the  fully  developed  sea  spectral  value. 
When  the  input  energy  and  the  dissipated  energy  balance  over  the 
whole  frequency  range,  the  sea  can  be  considered  to  be  a fully 
developed  sea.  There  is  some  question  as  to  whether  or  not  a 
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"fully  developed  sea"  exists  in  the  actual  ocean,  but  many  observed 
data  confirm  that  there  is  a stage  that  can  be  said  to  be  a fully  developed 
sea  and  that  the  spectral  components  are  proportional  to  the  f law 
in  the  high  frequency  region. 

When  the  growth  of  a spectral  component  is  considered,  this  energy 
dissipation  might  be  assumed  to  be  a function  of  the  ratio  of  the  spectrum 
at  a particular  moment  to  the  fully  developed  sea  spectrum  , and  as  a 
good  approximation,  equation  (2.  1)  may  be  expressed  by 


dt 


■S(f,  t,  x ) = [A(f,  u)  + B(f,  u ) S(f,  t,  jT)] 


[A(f,  u)  + B(f,ujS(f,t,x)] 


S(f,  t,  x ) 

S (f) 

00 
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(6.2) 


where  (f)  is  a fully  developed  sea  spectrum  for  that  wind  speed. 

In  equation  (6.2),  the  first  term  represents  the  energy  input  from  the 
atmosphere  to  the  waves  in  terms  of  a resonance  mechanism  and  an 
instability  mechanism,  and  the  second  term  includes  only  dissipation 
quantities  and  represents  the  energy  loss.  The  fully  developed  sea 
spectrum  form  proposed  by  Pierson  and  Moskowitz  (1964)  is  used. 

This  fully  developed  sea  spectrum  form  is  given  as 

,4 

(6.3) 
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where  a = 8.10  x 10  , |3  = 0.74,  and  wq  = g/u  . 

In  order  to  make  the  manipulation  of  (6.2)  easier,  another 


function  is  introduced  in  the  A term.  This  change  does  not  affect 
the  growth  very  much.  Then  (6.2)  is  changed  to  the  form 
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dt 


S(f,t,  3T)  = 
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The  solution  to  this  equation  for  a zero  initial  spectrum 
condition  and  an  infinite  fetch  is 


S(f,  t)  = ^sse^Liil  ^ + | *fegp  ■ (Bt a^iL , j2" 
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(6.5) 


where  A and  B are  represented  in  equations  (5.  1)  and  (5.  9) 

respectively.  If  t is  in  hours,  then  A and  B are  expressed  by 
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the  following  form  to  give  a spectrum  in  (m  sec). 
tt/2 
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If  the  sea  is  not  calm  initially,  equation  (6.8)  may  be  used. 


S(f,t)  = 


A[exp{B(t  + tQ)} 
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and  SQ  is  an  initial  spectral  component. 

The  spectral  growth  was  computed  for  wind  speeds  of  20, 

25,  30,  35,  40,  and  45  knots  by  using  (6.5).  The  results  are  shown  in 


figures  4(a)  through  (f)  as  a function  of  wind  duration  for  infinite 
fetch. 

If  the  wind  has  been  blowing  for  a long  time,  the  sea  is  then 
independent  of  time,  and  equation  (2.  1)  becomes 


^ S(f»  x)  = cg  S(f,  x)  = A(f,  u)  + B (f,  uj  S(f,  x)  (6 

From  x = c t,  the  spectral  growth  as  a function  of  fetch  was  also 
S 

computed  for  constant  wind  speeds,  and  the  results  are  shown 
in  figures  5(a)  through  (f). 

The  initial  condition  of  the  sea  was  zero  for  spectral  growths 
both  as  functions  of  fetch  and  duration.  From  this  initial  condition, 
the  friction  velocity,  u„. , cannot  be  obtained  from  equation  (5.6). 
Therefore,  the  initial  friction  velocity  u.,,  is  approximated.  The 
spectrum  for  the  first  time  step  is  computed,  and  then  u.,  for  the 

rv  — 

following  time  step  is  computed  from  the  values  of  c and  Hi  ob- 

3 

tained  from  the  spectrum  for  the  first  time  step. 

— /v 

As  the  waves  grow,  the  quantity  (u  - c)  influences  the 
friction  velocity  more  and  more.  That  is,  the  spectra  grow  in  the 

r*j  — 

low  frequency  band,  and  c approaches  the  wind  speed  u . Accord- 
ing  to  equation  (5.6),  u„.  decreases  as  (u  - c ) becomes  small. 
Therefore,  the  friction  velocity  has  a tendency  to  decrease  as  the 
waves  grow.  When  we  consider  the  following  form. 


= J c10  u 


where  c^q  is  a drag  coefficient  at  10  m above  the  sea  surface. 
c^q  also  decreases  as  the  waves  grow.  Kitaigorodskii  and  Volkov's 
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Figures  5c  and  5d  Growth  of  the 
spectrum  for  30  and  35  knot  winds 
with  fetch  in  nautical  miles. 
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Figures  5e  and  5f  Growth  of  the 
spectrum  for  40  and  45  knot  winds 
with  fetch  in  nautical  miles. 
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equation  was  verified  by  using  observations  representing  a rather 
wide  range  of  sea  conditions.  However,  none  of  these  observations 
(including  any  of  the  other  published  observations  of  the  drag  coeffi- 
cient Cjq  for  high  wind  speeds)  are  representative  of  sea  states 
for  fully  developed  seas.  When  equation  (5.6)  is  applied  to  spectra 
growing  in  these  final  stages,  c is  very  close  to  the  wind  speed, 
u,  and  the  corresponding  drag  coefficient  approaches  very  low  values. 
To  check  this  tendency,  a reasonable  lower  limit  of  0.0007  was  set 
at  the  ratio  u^/u  = ,s/c j q - Figures  4 and  5 were  computed  with 
this  assumption. 

In  these  figures,  the  spectra  develop  rapidly  over  the  high 
frequency  region  in  their  early  stages.  This  trend  happens  for 
all  wind  speed  cases.  For  frequencies  near  the  peak  of  the  fully 
developed  sea  spectra,  the  pace  of  the  spectral  growth  decreases. 
When  the  average  phase  speed  c comes  close  to  the  wind  speed 
u , the  friction  velocity  decreases,  and  the  growth  rate  of  equation 
(5.9)  also  decreases  in  a general  sense.  Therefore,  less  energy 
is  going  to  be  fed  in  from  the  atmosphere  to  the  waves  in  com- 
parison with  the  energy  input  in  the  earlier  stages.  In  this  sense, 
then,  the  growth  of  the  spectrum  depends  on  the  past  stages  of 
the  growth  and  on  the  frequencies  already  present  in  the  spectrum. 
This  can  be  interpreted  in  terms  of  the  fact  that  the  instability  mech- 
anism growth  is  not  so  dominant  in  the  low  frequency  region,  and 
therefore,  it  takes  a long  time  to  reach  the  fully  developed  stage. 

This  tendency  occurs  also  for  fetch  spectral  growth  with  respect  to 
distance. 


Sill!. 


The  spectral  form  of  the  partially  developed  sea  shows  a 
leading  forward  slope  of  the  spectra  that  is  not  very  steep  for  the 
early  stages  or  for  short  fetches.  Then  the  forward  face  becomes 
steeper  while  the  spectra  grow  in  the  lower  frequency  range.  As 
mentioned  previously,  this  is  the  effect  of  the  low  growth  rate  of  the 
instability  mechanism  at  those  frequencies.  In  figure  6,  the  spectra 
computed  from  the  theoretical  instability  growth  rate  (4.18)  are 
shown. 


7.  Comparison  with  other  proposed  spectra  and  spectral  growth  rates 

Many  theories  have  appeared  since  the  beginning  of  the  study  of 
actual  ocean  wave  generation.  First,  the  parameters  of  the  sea  were 
computed  by  the  non-dimensional  method  (Sverdrup  and  Munk,  1947; 
Bretschneider,  1957).  After  the  concept  of  the  spectra  was  introduced, 
wave  observations  were  made  from  the  spectral  point  of  view,  and 
many  spectral  forms  were  proposed  (Neumann,  1953;  Pierson, 
Neumann,  and  James,  1955;  Darbyshire,  1955,  1959;  Gelci,  Cazale, 
and  Vassal,  1957;  Bretschneider,  1959).  Directional  spectra  were 
observed  by  Cot£  et  al  (I960),  and  Longuet-Higgins  et  al  (1963). 

Walden  (1963)  gave  a review  of  the  differences  among  these  pro- 
posed wave  spectra.  Later  Pierson  and  Moskowitz  (1964)  proposed 
a form  for  fully  developed  sea  spectra  based  on  a similarity  theory. 

As  the  first  step,  fully  developed  sea  spectra  were  studied 
more  than  the  partially  developed  sea  spectra.  When  forms  of  the 
spectra  are  compared,  the  spectra  proposed  by  Darbyshire  and 
Gelci  et  al  (DSA  II)  show  smaller  spectral  densities  than  the  other 
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proposed  spectra  over  the  entire  range  of  frequencies.  It  is  not 
easy  to  compare  these  spectra  because  many  different  factors  such 
as  geographical  conditions  and  differences  in  instrumentation  are 
involved.  One  reason  for  the  discrepancies,  investigated  by  Pierson 
(1964),  was  the  wind  speed  variation  as  a function  of  height.  Pierson 
obtained  closer  agreement  among  the  different  proposals  by  Neumann, 
Wilson,  and  Pierson  and  Moskowitz. 

In  figure  7,  the  spectral  growth  for  a 40  knot  wind  at  19.5  m 
by  the  PNJ  method  is  shown.  Pierson  (1964)  compared  the  fully 
developed  sea  spectra  given  by  Neumann  (1953)  and  by  Pierson  and 
Moskowitz  (1964),  and  found  them  to  be  similar.  Partially  developed 
sea  spectra  in  this  study  grow  very  fast  for  0 to  15  hour  durations 
over  the  higher  frequency  range,  and  the  spectra  pass  the  peak  of 
the  fully  developed  sea  spectra  after  24  hours  from  a zero  initial 
spectrum.  According  to  the  PNJ  method  (H.  O Pub.  603),  the 
location  of  the  forward  edge  of  the  spectrum  is  shown  approximately, 
so  the  partially  developed  sea  spectra  in  figure  7 are  approximated 
by  their  leading  edges.  Generally,  the  partially  developed  sea 
spectrum  of  the  PNJ  method  develops  steadily  toward  the  lower 
frequency  range,  and  it  fills  the  fully  developed  sea  spectrum  after 
a 36  hour  duration.  The  spectrum  of  this  study  grows  fast  in  the 
early  stage,  but  does  not  attain  the  fully  developed  sea  spectrum 
easily  in  the  low  frequency  range.  If  the  final  stage  of  the  develop- 
ment is  neglected,  it  is  considered  that  the  spectrum  reaches  a fully 
developed  sea  after  40  hours. 

When  there  is  a big  difference  in  the  spectral  growth  theories. 
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the  growth  of  some  sea  state  parameters  must  be  investigated.  The 
significant  wave  height  Hi  as  a function  of  the  total  energy  of  the  sea 

3 

maybe  an  appropriate  one  to  check  wave  growth.  The  significant 
wave  height  is  plotted  as  a function  of  duration  or  fetch  in  figures 
8(a)  and  8(b).  This  significant  wave  height  = 2.83'\/E  , where  E 

3 

is  twice  the  variance  of  the  spectrum.  In  the  same  figures,  the  growth 

of  Hx  obtained  by  Sverdrup  and  Munk  (1947)  and  PNJ  (1955)  are 

3 * 

also  shown  for  wind  speeds  of  20,  30,  and  40  knots  at  19.5  m.  The 
growth  of  the  significant  wave  height  proposed  by  Sverdrup  and  Munk 
fits  very  well  for  20  and  30  knot  winds  in  the  duration  growth  graph 
8(a).  The  wind  speed  applied  was  corrected  to  the  value  at  the 
height  of  19.5  m by  using  the  drag  coefficient  proposed  by  Deacon 
and  Webb  (1962)  from  10  m for  Sverdrup  and  Munk1  s growth,  and 
from  7.5  m for  the  PNJ  method  growth.  However,  the  fetchwise 
growth  of  this  scudy  is  more  like  the  PNJ  method  than  the  Sverdrup- 
Munk  method. 

The  growth  rate  of  Hx  by  the  PNJ  method  increases  more 

3 

slowly  than  this  proposed  growth.  This  is  clear  from  the  partially 
developed  sea  growth  in  figures  4 and  5.  Similar  figures  for  40 
knot  winds  significant  wave  height  growth  are  shown  by  Walden 
(1963)  and  Wilson  (1965)  with  many  other  proposed  growth  curves. 
However,  the  height  at  which  wind  speed  was  measured  is  not 
clearly  specified  in  these  growth  rates.  There  are  many  variations 
in  these  figures.  Generally,  there  are  two  growth  rate  groups.  One 

*Early  definitions  of  significant  height  were,  of  course,  not  this 
precise. 
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group  shows  a fast  growing  sea,  and  the  other  shows  relatively  slow 
growth. 

One  of  the  most  distinctive  points  is  the  growth  of  the  significant 
wave  height  proposed  by  Darbyshire.  According  to  his  study,  the  sea 
reaches  the  fully  developed  state  easily  for  any  wind  speed.  100  to 
150  nautical  miles  as  a fetch,  and  10  to  15  hours  as  a duration  are 
said  to  be  enough.  On  the  contrary,  500  or  more  nautical  miles  and 
36  hours  are  needed  for  the  40  knot  wind  speed  spectra  with  the  PNJ 
method. 

There  may  be  two  reasons  for  considering  a short  fetch  re- 
quirement for  fully  developed  seas.  One  reason  is  as  follows.  When 
a wind  field  itself  is  moving  with  the  same  velocity  as  the  propagation 
speed  of  the  wave  energy  in  the  direction  of  the  waves,  then  the  waves 
can  grow  considerably  higher  even  for  a short  apparent  fetch.  The 
second  reason  is  the  geographical  location  of  the  observation  point.  If 
the  observation  point  is  exposed  to  waves  that  are  propagating  from 
other  areas,  it  is  hardly  possible  for  that  location  to  be  glassy  calm. 

The  weather  ship  stations  in  the  North  Atlantic  can  be  considered  to 
be  locations  such  as  this.  While  the  sea  always  has  some  energy,  an 
initial  excitation  due  to  the  resonance  mechanism  is  no  longer  needed, 
and  so  the  sea  grows  in  an  exponential  manner  immediately  and  the 
duration  required  to  reach  a fully  developed  state  becomes  small. 

In  an  actual  ocean,  the  two  reasons  mentioned  above  may  occur 
simultaneously. 

Suppose  there  exists  some  background  sea.  The  spectrum 
grows  faster  and  so  does  the  wave  height.  To  investigate  this  back- 
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ground  sea  effect,  the  background  sea  level  was  assumed  white 
and  given  a variance  of  a 20  knot  wind  sea.  Furthermore,  before 
strong  winds  of  40  or  50  knots  blow,  some  weaker  wind  has  already- 
occurred,  and  a spectrum  that  is  equivalent  to  the  fully  developed 
sea  spectrum  for  10  knots  is  supposed  to  exist  in  addition  to  the 
white  noise.  The  larger  of  the  values  from  either  the  white  noise  or 
the  10  knot  wind  fully  developed  sea  was  used.  The  white  noise  vari- 

2 _ i 

ance  was  0.011  (ft)  for  a frequency  band  of  1/180  sec  . The 
growths  of  the  partially  developed  sea  spectra  for  30  knots  and 
40  knots  were  computed  as  shown  in  figure  9(a)  and  figure  9(b)  under 
the  background  sea  assumption.  The  background  level  used  is  also 
indicated  in  figure  9(a),  but  not  in  9(b),  because  the  level  is  so  small 
that  it  cannot  be  seen.  Naturally,  the  spectrum  started  from  the 
background  sea  grows  faster  and  it  reaches  a fully  developed  state 
easily  because  the  growth  in  the  low  frequencies  is  quite  rapid. 

The  discussion  of  the  growth  for  a very  short  fetch  or  very 
short  duration  may  have  no  meaning  in  this  study.  Even  if  the 
spectral  growth  for  a particular  component  is  the  same,  the  different 
ones  proposed  have  different  spectral  forms.  So  the  growth  behavior 
of  the  significant  wave  height  is  also  different.  However,  the  above 
reasons  for  a faster  growth  in  some  theories  may  be  valid. 

One  of  the  points  that  should  be  noted  in  this  model  is  that  a 
particular  spectral  component  shows  a different  growth  behavior  even 
for  the  same  significant  wave  height  and  the  same  wind  speed.  The 
spectral  growth  depends  on  the  friction  velocity  and  not  simply  on  the 
wind  speed,  and  so  the  spectrum  grows  at  different  rates,  even  for 
the  same  wind  and  the  same  wave  height. 
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For  example,  consider  the  two  different  spectra  described 
below.  Suppose  one  spectrum  has  a steep  forward  slope  that  covers 
quite  a narrow  frequency  band  width.  The  other  spectrum  has  a 
gentle  forward  slope  that  covers  a wider  frequency  band.  Then 
there  exist  two  spectral  components  whose  values  are  the  same. 
However,  the  growths  of  these  two  spectral  components  are  different 
because  there  are  different  friction  velocities  according  to 
equation  (5.6). 


8.  Application  to  forecasting 


One  of  the  purposes  of  ocean  wave  studies  is  to  help  forecast 
waves,  and  one  object  of  this  study  is  to  imp.j/e  wave  forecasting 
methods.  Therefore,  computations  by  using  this  model  must  be 
compared  with  sea  conditions  that  actually  existed.  For  that  purpose, 
it  is  desirable  to  hindcast  the  waves  over  a wide  area  and  pick  one 
particular  geographical  point  to  check  the  computed  result. 

Since  the  high  speed  computer  has  become  available,  the  study 
of  wave  forecasting  has  been  possible  over  wide  ocean  areas,  and 
has  shown  remarkable  progress.  Baer  (1962)  studied  wave  fore- 
casting methods  in  the  North  Atlantic,  and  used  519  grid  points 
which  were  roughly  120  nautical  miles  apart  over  the  whole  North 
Atlantic  Ocean.  At  that  time,  the  change  of  the  quantity,  E,  was 
tabulated  per  unit  time  interval  by  using  the  PNJ  method,  and  the 
growth  of  the  following  time  steps  was  obtained  from  that  growth 
table.  Many  parts  of  the  program  were  revised  later,  and  the  de- 
tails of  the  numerical  procedure  were  given  in  other  publications 
(Pierson,  1964;  Pierson,  Tick  and  Baer,  1966;  Bunting,  1966). 


The  essence  of  this  process  is  described  below. 

First  of  all,  the  wind  field  must  be  determined.  The  regression 
equation  given  by  Thomas  ell  and  Welsh  (1963)  is  used  to  guess  the 
first  wind  field  from  the  surface  pressure  field.  This  estimate,  then, 
is  corrected  by  using  the  observed  wind  data  obtained  by  weather  ships, 
navy  ships,  and  other  commercial  ships  with  priority  in  this  order. 

The  wind  speeds  are  referred  to  a 19.5  m height  which  is  the  British 
weather  ships'  anemometer  height  by  assuming  the  wind  profile  to  be 
logarithmic.  The  wind  profile  is  obtained  in  the  JNWP  grid  system, 
and  is  interpolated  into  the  519  grid  point  system  mentioned  above 
(Baer,  1962). 

The  growth  of  the  spectra  for  three -hour  intervals  is  computed 
under  the  wind  field  obtained.  Equation  (6.3)  is  used  for  fully  developed 
sea  spectra.  When  the  spectra  at  each  time  are  below  the  fully  developed 
spectral  component  value,  further  growth  is  allowed.  The  spectral 
growth  model  of  this  study  is  applied  to  a growth  instead  of  the  E 
growth  table.  The  forecasts  (or  hir.dcasts  in  this  study)  of  the  spectra 
are  made  for  fifteen  frequency  bands  to  make  computation  faster. 

The  following  frequency  bands  were  chosen  [unit  is  cycles  per  second]. 


0.036  - 0.042, 


0.042  - 0.047  , 


0.047  - 0.053  , 


0.053  - 0.G58, 


0.058  - 0.064  , 


0.064  - 0.069  » 


0.069  - 0.075, 


0.075  - 0.080  , 


0.080  - 0.086  , 


0.086  - 0.097  , 


0.097  - 0.108  , 


0.108  - 0.125  , 


0.125  - 0.142  , 


0.142  - 0.164  , 


0. 1 64  - oo 
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In  the  fifteenth  band,  the  sea  is  estimated  to  be  always  fully 
developed  since  the  high  frequencies  easily  reach  the  fully  developed 
sea  spectrum.  Thus,  the  fully  developed  spectral  value  is  applied 
to  that  band  instead  of  computing  a growth.  The  functions  for  A in 
equation  (6.6)  and  for  B in  equation  (6.7),  must  be  obtained  first, 
and  the  A term  is  simply  a function  of  wind  speed  and  frequency. 

The  B function  depends  on  a friction  velocity  u..  , and  u_  should 
be  computed  from  equation  (5.6).  Thus,  before  computing  the 
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spectral  growth,  the  sea  conditions,  and  c , must  be  obtained 

3 

/v 

for  the  grid  point  concerned.  When  an  average  c is  computed,  the 
2 

quantity  f S(f)  is  needed.  The  center  frequency  f^  of  the  band  is 
2 

used  and  fcS(fc)  is  applied  even  for  the  broader  frequency  band 
except  the  fifteenth  frequency  band  (0.164  cps  ~ °o  ). 

The  friction  velocity,  , for  each  grid  point  is  obtained 
from  the  initial  sea  state,  and  the  spectral  component  growth  for 
three  hours  is  computed  by  using  this  model.  The  spectrum  ob- 
tained in  such  a way  is  a one -dimensional  spectrum,  so  the  equation 
derived  by  the  SWOP  project  (Cot£  et  al,  I960)  is  used  to  obtain  the 
directional  spectra.  The  formula  is 

4 

F(w,  0,  u)  = - [1  + (0.50  + 0.82e"^wu/g)  ) cos  20 


t0.32e'^uu/g^  cos  40] 


for  - ^ < 0 < , and  F(u,  0,  u)  = 0,  elsewhere.  0 is  the  angle 

between  the  wind  direction  and  the  direction  of  the  wave  component. 
The  directional  spectra  are  computed  for  30°  angular  bands.  As 
a result,  the  spectra  are  shown  by  180  spectral  components 
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(15  frequency  bands)  x (12  directional  bands).  The  spectral  values  are 
propagated  at  group  velocity  in  their  appropriate  direction. 

While  the  wave  grows  and  propagates,  the  wave  may 
also  dissipate  in  its  travels.  The  spectral  component  can  be  con- 
sidered not  to  receive  any  energy  but  to  dissipate  when  the  wave 
moves  against  the  wind.  In  this  forecasting  procedure,  the  dissipation 
rate  is  assumed  to  be  a function  of  the  total  energy  of  the  wind  sea 
traveling  within  90°  to  the  wind  direction,  and  the  fourth  power  of  its 
component's  frequency. 
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where 


Sjj(T,  0^)  = spectral  component  after  dissipation 
SQ(f.,  0^)  = spectral  component  before  dissipation 
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Thus,  the  spectral  components  computed  at  a grid  point  transfer  some 
energy  to  the  other  grid  points,  receive  energy  from  the  other  grid 
points,  and  combine  to  form  a new  sea  state.  The  next  time  step 
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computation  is  repeated  from  this  new  sea  condition,  and  the  result  is 
recorded  every  six  hours. 

Wave  hindcasts  for  December  1959  were  made  for  the  North 
Atlantic  as  a test.  In  the  second  half  of  December  1959,  three  big 
storms  attacked  the  northern  part  of  that  ocean.  The  meteorological 
and  sea  conditions  are  discussed  in  detail  by  Bretschneider  et  al  (1962). 
Within  half  a month,  significant  wave  heights  of  forty  feet  were  ob- 
served twice,  and  a wind  speed  of  more  than  60  knots  was  observed 
at  times  on  the  weather  ship  located  at  "J".  The  corresponding 
changes  of  the  wave  height,  which  accompanied  the  changes  of  the 
wind  speed,  can  be  thought  to  be  the  best  sample  for  checking  this 
growth  model  by  using  the  forecasting  procedure  mentioned  above. 

The  wind  field  used  in  this  test  was  derived  by  the  regression 
equation  from  the  pressure  pattern,  and  corrected  by  the  ship's  data. 
Thus  there  were  some  differences  between  the  actual  wind  field  and 
the  input  wind  field,  and  also  the  grid  point  chosen  was  not  at  the 
same  location  as  the  weather  ship.  The  observed  wind  velocities 
and  the  nindcasted  wind  velocities  are  listed  in  Appendix  III.  Atmo- 
spheric stability  was  not  taken  into  account  in  this  study,  so  the  air- 
sea  temperature  difference  effect  was  neglected. 

The  computation  starts  from  a zero  initial  sea  condition,  but 
there  are  enough  days  to  produce  a sea  state  that  is  close  to  the 
actual  sea  condition  before  16  December.  The  computation  began 
eight  days  prior  to  16  December.  The  time  history  of  the  hindcasted 
significant  wave  height  at  the  two  grid  points  closest  to  weather  ship 
J is  shown  in  figure  10.  The  relative  position  to  J is  seen  in  that 
figure. 
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Although  there  are  some  differences  between  the  observed 
and  hindcasted  significant  wave  heights,  these  hindcasted  values  follow 
the  observed  ones  fairly  well.  The  bias  of  the  wave  height  is  + 0.4  m 
(+1.3  ft)  and  the  RMS  error  is  1.6  m (4.9  ft)  for  grid  point  72,  and 
+0.7  m (2.1  ft)  and  1.4  m (4.3  ft),  respectively,  for  grid  point  73. 

The  purpose  of  this  wave  forecasting  study  is  not  only  to  pre- 
dict significant  wave  height,  but  also  to  predict  spectra  close  to  the 
actual  ones  observed.  Therefore,  the  investigation  must  be  extended 
to  the  form  of  the  spectra  or  to  the  spectral  components  changing 
with  time.  The  study  to  check  >w  the  ocean  wave  spectra  can  be 
reproduced  is  important  for  practical  purposes  such  as  ship 
motions  in  waves.  The  one-dimensional  spectra  were  obtained  to 
study  the  growth  and  decay  of  the  spectral  components. 

Although  fifteen  frequency  bands  were  computed,  only  six 
low  frequency  bands  are  shown  in  figure  11.  The  spectral  com- 
ponents in  the  higher  frequency  bands  showed  growth  and  decay  that 
agreed  well  with  the  observations,  and  a problem  remained  in  the 
low  frequency  range  only.  When  the  old  growth  model  (Inoue,  1966) 
was  used,  there  was  no  reproduction  of  spectral  components  that 
appeared  in  the  low  frequencies.  Bue  in  this  model,  even  in  the 
lowest  frequency,  f = 0.039  (cps),  the  spectral  component  was 
reproduced  quite  well.  This  low  frequency  component  prediction 
is  considered  to  be  the  effect  of  the  undulatory  turbulent  flow  induced 
mainly  by  the  waves.  The  hindcasted  components  seem  to  grow  a 
little  too  fast,  although  they  do  not  decrease  fast  enough  after  the 
storm  for  f = 0.055,  0.050,  and  0.045  cps.  The  results  could  not 
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be  further  refined  without  investigating  the  actual  wind  field  in 
detail  over  the  whole  North  Atlantic.  As  far  as  the  low  frequency 
part  of  the  spectrum  is  concerned,  the  growth  rate  of  this  study  im- 
proved the  forecasting  results. 

Hindcasts  of  ocean  waves  were  also  made  for  20  November 
through  30  November  1961.  During  this  period,  wave  observations 
were  made  at  the  Argus  Island  tower  and  analyzed  by  Pickett  (1962).  The 
meteorological  and  wave  situations  were  also  discussed  by  De  Leoni- 
bus  (1962)  as  well  as  the  wave  hindcasts  obtained  by  using  the  PNJ 
method.  The  hindcasted  significant  wave  heights  for  the  two  grid 
points  nearest  to  the  tower  are  shown  in  figure  12  along  with  the 
observed  significant  heights.  The  procedure  for  the  computation  is 
the  same  as  the  December  1959  case  except  for  the  wind  field  used, 
which  was  analyzed  manually  [provided  by  L.  Moskowitz*].  Al- 
though there  are  differences  at  the  two  grid  points  due  to  the  com- 
plicated weather  situation  for  the  first  half  of  the  period,  the  most 
interesting  hindcasts  occurred  on  25,  27,  and  28  November.  These 
rapid  wave  height  changes  have  never  been  reproduced  before.  The 
bias  and  the  RMS  error  for  grid  point  229  are  respectively, 

-0.53  m (-1.59  ft)  and  1.25  m (3.76  ft);  for  grid  point  230  they  are 
-0.17  m (-0.57  ft)  and  0.95  m (3.10  ft).  Generally  speaking,  the  hind- 
casts from  this  study  are  more  sensitive  to  the  meteorological  and 
sea  conditions  in  comparison  to  the  growth  table  method  and  the  old 
model  (Inoue,  1966).  This  is  one  of  the  distinguishing  characteristics 
of  this  procedure. 


*U.  S.  Naval  Oceanographic  Office. 
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9.  Conclusions 


Observations  showed  an  irregular  trend  in  the  figure  of 
instability  mechanism  growth,  as  shown  in  figure  3.  This  "bump" 
had  not  been  explained  sufficiently.  Although  there  are  discrepancies 
between  observed  growth  rates  and  theoretical  values  in  the  high  and 
also  low  ranges  of  u^/c  , the  new  theory  proposed  by  Phillips  can 
explain  the  existence  of  the  "bump"  in  figure  3.  The  growth  rate 
is  a function  of  u..{/c,  and  the  friction  velocity,  u_,_ , should  be  deter- 
mined. The  relation  between  the  roughness  parameter  and  the  sea 
state  was  found  by  Kitaigorodskii  and  Volkov,  and  this  relation  was 
applied  to  the  computation  of  the  friction  velocity. 

Therefore,  the  friction  velocity  is  a function  of  the  sea  state, 
or  the  average  phase  speed  and  the  wave  height,  and  so  is  the  growth 
rate,  B/f  . Even  for  the  same  wind  speed  at  19.5  meters,  the 
growth  rate  for  a particular  component  is  not  the  same,  but  depends 
on  the  sea  state.  Physically,  this  can  be  interpreted  to  mean  that 
less  energy  can  be  fed  in  by  an  instability  process,  while  the  waves 
are  growing  and  moving  with  the  wind. 

The  proposed  sea  spectra  growth  shows  faster  development 
in  an  early  growing  stage  (or  short  fetch),  and  slower  growth  in  the 
last  stages.  It  takes  from  35  to  40  hours  or  from  700  to  1000  nauti- 
cal miles  to  reach  a fully  developed  sea  for  a 40  knot  wind.  This 
growth  rate  of  the  significant  wave  height  is  almost  the  same  as  the 
results  of  Sverarup  and  Munk  for  duration,  but  for  fetch  the  growth 
is  slower.  The  durations  and  fetches  required  to  reach  fully  de- 
veloped seas  seem  larger  in  comparison  to  the  other  proposed  values 


such  as  those  of  Darbyshire.  The  discrepancies  maybe  explained 
by  the  existence  of  a background  sea. 

The  formulas  obtained  in  this  study  were  used  in  a wave 
forecasting  computer  program.  The  hindcasted  significant  wave 
height  shows  good  agreement  with  the  observed  heights  at  weather 
station  J in  the  North  Atlantic  for  the  second  half  of  December  1959. 

Spectral  components  were  also  investigated  for  low  frequency 
ranges  which  had  not  been  properly  described  in  an  earlier  model. 

The  spectral  components  in  that  region  also  showed  better  growth. 
Although  this  is  a linear  model  basically  and  the  atmospheric 
stability  was  assumed  to  be  neutral,  the  results  can  be  said  to  be 
satisfactory. 
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Appendix  II.  Air-sea  temperature  difference  (°C) 
British  weather  ships  (Moskowitz,  Pierson  and  Mehr) 


JHC  38-39 

Apr.  8,  1955 

(-0.5)  ~ (-0.4) 

JHC  73-74 

May  11,  1956 

(-1.5)  ~ (-0.2) 

DL  37-38 

Oct.  28,  1958 

(-0.9)  ~ (-3.9) 

JH  10-11 

Mar.  29,  1959 

(-1.1)  ~(-2.1) 

DL  55-56 

Nov.  9,  1959 

(-5.8)  ~ {-7.0} 

JHA  9-10 

Dec.  17,  1959 

(-3.8)  ~ (-3.1) 

DL  77-78 

Jan.  30,  I960 

(-0.6)  ~ (-1.7) 

Airborne  altimeter  (Barnett  and  Wilkerson,  1966} 

(-6C)  ~ (-10°) 

Towed  raft  (Snyder  and  Cox,  1966) 

Almost  neutral  condition. 

Largest  difference  (-0.7) 
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Appendix  III.  Wind  velocity  for  the  test  run  period, 
December  1959. 


Weather 

ship 

Grid  point  72 

Grid  point  73 

Date 

time 

(GMT) 

Speed 

(kts) 

Dir. 

(°) 

Speed 

(kts) 

Dir. 

(°) 

Speed 

(kts) 

Dir. 

(°) 

Dec.  16 

00 

10 

340 

25.3 

270 

25.7 

278 

06 

04 

280 

21.1 

315 

22.1 

246 

12 

23 

170 

28.3 

209 

28.0 

194 

18 

32 

250 

37.9 

243 

36.9 

234 

Dec.  1 V 

00 

60 

250 

34.7 

270 

35.8 

257 

06 

48 

280 

41.9 

268 

44.8 

252 

12 

48 

290 

48.0 

286 

46.8 

276 

18 

48 

290 

48.7 

293 

50.4 

283 

Dec.  18 

00 

40 

320 

21.8 

297 

33.6 

295 

06 

20 

280 

16.3 

277 

30.7 

288 

12 

26 

180 

31.2 

160 

18.8 

171 

18 

26 

210 

34.2 

171 

32.5 

165 

Dec.  19 

00 

30 

220 

22.5 

216 

23.1 

207 

06 

18 

250 

27.6 

221 

28.8 

215 

12 

28 

240 

26.0 

240 

30.0 

232 

18 

17 

250 

30.7 

233 

31.6 

224 

Dec.  20 

00 

28 

270 

24.2 

264 

27.7 

252 

06 

36 

280 

24.3 

260 

29.7 

253 

12 

28 

260 

29.9 

255 

31.0 

259 

18 

16 

250 

27.1 

262 

30.7 

251 

Dec.  21 

00 

28 

270 

26.8 

272 

28.3 

258 

06 

27 

310 

16.5 

269 

21.8 

259 

12 

07 

270 

7.9 

239 

15.7 

264 

18 

23 

170 

18.2 

178 

14.0 

187 

Weather  ship  Grid  point  72  Grid  point  73 


Date  time 
(GMT) 

Speed 

(kts) 

Dir. 

<°) 

Speed 
(kt  s) 

Dir. 

Speed 

(kts) 

Dir. 

(°) 

Dec.  22  00 

35 

220 

32.8 

220 

27.9 

201 

06 

40 

230 

47.2 

237 

43.1 

231 

12 

40 

250 

37.8 

251 

37.4 

230 

18 

50 

270 

38.2 

244 

39.2 

224 

Dec.  23  00 

55 

270 

45.5 

269 

45.9 

257 

06 

43 

280 

45.2 

257 

44.1 

247 

12 

45 

290 

41.  7 

287 

40.4 

269 

18 

(21Z) 

40 

300 

41.0 

292 

41.0 

281 

Dec.  24  00 

34.8 

293 

38.6 

293 

06 

32 

290 

22.9 

285 

32.7 

287 

12 

20 

270 

7.5 

249 

20.0 

275 

18 

18 

160 

12.1 

144 

6.3 

216 

Dec.  25  00 

43 

130 

24.7 

129 

32.8 

128 

06 

21 

160 

14.3 

166 

24.3 

J 59 

12 

19 

230 

9.8 

244 

17.0 

206 

18 

25 

210 

17.2 

264 

26.1 

224 

Dec.  26  00 

37 

250 

22.2 

273 

28.9 

241 

06 

23 

270 

26.8 

287 

31.4 

259 

lc 

22 

290 

21.1 

293 

24.9 

259 

18 

21 Z 

18 

360 

24.5 

295 

27.4 

271 

Dec.  27  00 

03 

23 

360 

19.3 

295 

22.6 

274 

06 

09 

28 

320 

23.3 

298 

26.4 

293 

12 

15 

25 

330 

26.6 

299 

29.8 

303 

18 

21 

23 

300 

20.7 

289 

27.2 

285 
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